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Quantum spin clusters with dominant antiferromagnetic Heisenberg exchange interactions typi- 
cally exhibit a sequence of field-induced level crossings in the ground state as function of magnetic 
field. For fields near a level crossing, the cluster can be approximated by a two-level Hamiltonian at 
low temperatures. Perturbations, such as magnetic anisotropy or spin-phonon coupling, sensitively 
affect the behavior at the level-crossing points. The general two-level Hamiltonian of the spin sys- 
tem is derived in first-order perturbation theory, and the thermodynamic functions magnetization, 
magnetic torque, and magnetic specific heat are calculated. Then a magneto-elastic coupling is in- 
troduced and the effective two-level Hamilitonian for the spin-lattice system derived in the adiabatic 
approximation of the phonons. At the level crossings the system becomes unconditionally unstable 
against lattice distortions due to the effects of magnetic anisotropy. The resultant magneto-elastic 
instabilities at the level crossings are discussed, as well as the magnetic behavior. 
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I. INTRODUCTION AND CONCLUSIONS 

The issue of level crossings (LCs) in quantum spin 
clusters is of broad interest. The most elementary ex- 
ample is probably that of an effective spin-1/2 with an 
avoided LC in a magnetic field, where the tunneling prob- 
ability at constant field-sweep rate is described by the 
Landau-Zener-Stueckelberg model, 1 Experimental real- 
izations are, for instance, the single-molecule magnets 
Mni2 or Fes^ 

In antiferromagnetic (AFM) spin clusters, LCs also 
emerge from the interplay of Heisenberg exchange inter- 
actions in the cluster and the Zeeman term, which leads 
to a sequence of LCs in the ground state as function of 
magnetic field ,2iii§iSi2i2iS More precisely, spin clusters are 
considered which are modeled by the microscopic spin 
Hamiltonian 



H 



N 



Sj + /j, B gS ■ B + Hi 



(1) 



consisting of the Heisenberg exchange term and the Zee- 
man term (N is the number of spin centers in the clus- 
ter and S = 'J2 i Si is the total spin operator). The 
term Hi contains all remaining relevant terms, such as 
magnetic anisotropy (single-ion anisotropy, anisotropic 
exchange, dipole-dipole interactions, etc.) or couplings 
to the environment (phonons, intercluster magnetic cou- 
plings, nuclear spins, etc.), but should be small in the 
sense of a perturbation. Then the eigenstates of H can 
be classified by the total spin quantum numbers S and 
M. The AFM Heisenberg interactions typically lead to 
an energy level scheme, where the energy of the lowest 
state for each value of S > Sq increases according to 
E(So) < E(S + 1)< ...< E{S max ) {So is the ground- 
state spin, and S ma x — Y^l As function of field, the 
Zeeman splitting produces a sequence of LCs, at which 



the ground state changes from \S = Sq,M = —So) to 
\S = S + 1,M = -So-1), to \S = S + 2,M = -S -2), 
and so on4 For magnetic fields close to the LC condi- 
tions, the system becomes almost degenerate, and is well 
described as a two-level system at low temperatures. 

Studies on this type of LCs are interesting from sev- 
eral perspectives, and numerous experimental examples 
are known both in the area of molecular nanomagnets 
and inorganic magnetic compounds. For instance, the 
determination of the LC fields from field-dependent mea- 
surement of thermodynamic quantities, such as the mag- 
netization or magnetic torque, directly yields the ener- 
gies E(S + 1) - E(S ), E(S + 2) - E(S ), etc., of the 
next higher spin levels. This "thermodynamic" spec- 
troscopy often enables a precise determination of the 
exchange coupling constants in the clusteri2i£i5i£iZi2i2 In 
some molecular nanomagnets, such as the AFM wheels, 
the energy of the lowest spin levels increases quadratically 
with S according to E(S) oc S{S + l), 3 ! 10 ' 11 ' 12 which re- 
flects the quantized rotation of the Neel vector in these 
clusters^ii^^ 

Interesting phenomena also emerge from the fact, that 
near a LC the system is rather susceptible to the pertur- 
bations represented by H±. Among the molecular nano- 
magnets, the AFM wheels 3 -^ and related systems, such 
as the Mil- [3 x 3] grid or the modified wheels jiSiliaJiSi 
have become prototypical examples in this regard. The 
most important term in Hi usually is the single-ion 
anisotropy. In the wheels, however, because of their 
nominally high symmetry, the LCs would remain true 
LCs despite this perturbation ; 22 : 23 but with small devi- 
ations from ideal cyclic symmetry spin terms such as a 
Dzyaloshinsky-Moriya interaction may mix the two levels 
at a LC, thus leading to avoided LCs. 23,24 Indications of 
this has been reported for some ferric wheels i 23 ' 24 : 25 In 
a priory less symmetric molecules the symmetry restric- 
tions are absent, and the single-ion anisotropy becomes 
efficient in mixing the spin levels at the LCs. 19 : 26 : 27 : 28 
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The result can be dramatic. For instance, the magnetic 
torque at low temperatures may display a pronounced os- 
cillatory field dependence, i.e., show quantum magneto- 
oscillations. This phenomenon has been observed in the 
Mn-[3 x 3] gridj^ and two modified wheels. 2 - Since lev- 
els with different values of S are mixed, the total spin S 
exhibits quantum fluctuations at the LCs, which could 
enable coherent oscillations of the total spin if dissipa- 
tion is weak3i Moreover, the torque oscillations observed 
in Mn-[3 x 3] provided the experimental demonstration 
of the so-called tower of states in a finite AFM square 
latticed 

H\ also may include perturbations such as a coupling 
to the nuclear spins or the phonon system. In fact, 
the LCs in the molecular wheels were intensively stud- 
ied by nuclear magnetic resonance (NMR) | 23 i 29 i 30 i 31 i 32 
Also effects of the spin-phonon coupling were observed 
in the wheels NaFe 6 and Fe^ ) 25 ' 33 and other clusters* 3 ^ 
where it resulted in magnetic butterfly-hysteresis at the 
LCs due to a phonon-bottleneck effect. From a general 
perspective, at the LCs the clusters represent two-level 
systems with dissipation. 35 They open attractive experi- 
mental opportunities for exploring this model, since, for 
instance, the strength of the mixing at the LCs and 
hence the influence of dissipation may be tuned by exter- 
nal parameters such as the orientation of the magnetic 
field^2L2£ 

Other interesting examples are found in inorganic 
compounds, for instance, in the materials characterized 
as weakly-interaction dimers. Prototypical systems are 
TlCuCl 3 and BaCuSiaOg^i 3 ! The dominant AFM in- 
teractions in the dimers lead to a S = dimer ground 
state and a higher lying S = 1 level, which is split in a 
magnetic field resulting in a LC, exactly as described be- 
fore. Weaker magnetic exchange interactions between the 
dimers, however, drastically modify the behavior of the 
system for fields close to the LC; they give rise to various 
field-induced magnetic long-range ordering phenomena, 
such as the Bose-Einstein condensation of triplets and 
dimensional reduction at quantum critical points! 38 ' 39 ! 40 

The (near) degeneracy of the levels at the LCs suggests 
also the possibility of spontaneous structural distortions 
of the cluster at the LC fields, driven by the interaction 
of the spin system with the phonons. In fact, in view of 
the resemblance of an AFM molecular wheel with an ID 
AFM quantum spin chain, a spin-Peierls type of effect 
seems to be obvious. However, general arguments imply 
that for finite spin clusters, which exhibit a gap in their 
excitation spectrum, the energies should vary quadrat- 
ically with the distortion coordinate. Hence, a sponta- 
neous structural instability would be conditional, and not 
likely to occur in real materials. Furthermore, a recent 
theoretical study on AFM Heisenberg rings demonstrated 
that, for Si > 1/2, a spontaneous distortion would be 
of first order and so strong, that it would disrupt the 
molecule4i 

Recently it has been shown, however, that the mag- 
netic anisotropy terms in Hi in fact give rise to an un- 



conditional instability of the clusters against spontaneous 
distortions at the LCs^ 2 - The argument is based on a phe- 
nomenological model, which describes the spin system by 
a two-level Hamiltonian and treats the phonons in the 
adiabatic limit (the effective Hamiltonian of the coupled 
spin-phonon system will be called Hsp). In the molecu- 
lar wheel CsFes, phase transitions at the LC points were 
recently observed^ 2 - Since the crystal structure as well 
as other arguments disapproved the presence of sufficient 
magnetic interactions between the clusters, the findings 
were interpreted as magneto-elastic instabilities at the 
LCs. The experimental data in fact were consistent with 
the predictions of the phenomenological model Hsp- 

The primary goal of this work is to present a coherent 
derivation of Hsp and the field-induced magneto-elastic 
instabilities, and to discus the resultant magnetic behav- 
ior (extending a previous brief report^ 2 -). As a prepara- 
tory step, the two-level Hamiltonian of the spin system is 
derived (Sees. II and III), and the thermodynamic mag- 
netization, magnetic torque, and magnetic specific heat 
calculated (Sec. IV). Then the magneto-elastic coupling 
is introduced and the two-level Hamiltonian Hsp estab- 
lished (Sec. V), and the magneto-elastic instabilities dis- 
cussed (Sec. VI). The subject of field-induced LCs has 
been investigated intensely in the last years by several 
authors. Accordingly, several pieces of the results were 
obtained before. 

The results should be useful in several regards. The 
thermodynamic results provide insight into the different 
magnetic behavior of, e.g., AFM wheels and the Mn- 
[3 x 3] grid (in the wheels, the sequence of LCs leads 
to staircase-like magnetization and torque curves; in the 
grid, it leads to oscillations in the torque). Furthermore, 
they provide a handy method to parameterize experi- 
mental data, and to extract from them key parameters, 
such as the level-mixing strength. The presented model 
of the magneto-elastic instability has all the drawbacks 
and virtues of a phenomenological model. It cannot be 
expected to provide a full quantitative description of a 
specific material (e.g., of CsFes), or to answer specific 
microscopic details, such as which distortion mode is in- 
volved. However, it is the most simple model which de- 
scribes the basic mechanism (it actually represents the 
mean- field model). As it is often with such models, they 
are very useful for a qualitative understanding of the phe- 
nomenon, and discussion of observed behavior. 

A connection between the weakly-interacting dimer 
compounds and AFM molecular wheels is finally noted. 
The low-lying excitations in AFM wheels are well de- 
scribed by a dimer model, where the two spins corre- 
spond to the total spins on each AFM sublattice i 12 ' 13 ' 15 
Hence, with magnetic interactions between the clus- 
ters, AFM wheels would perfectly mimic a weakly- 
interacting dimer system, and show the very same field- 
induced phenomena. In the AFM wheels, however, the 
magnetic intercluster interactions are usually negligible, 
which suggested the magneto-elastic origin of the ob- 
served field- induced phase transitions in CsFes ^ In the 
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weakly-interacting dimer compounds the observed field- 
induced phase transitions are undoubtly of magnetic ori- 
gin, but also clear indications of magneto-elastic cou- 
plings have been reported i 43 i 44 ' 45 In particular, NMR ex- 
periments showed that the field-induced magnetic order- 
ing in TICUCI3 is accompanied by a simultaneous lattice 
deformation^ 3 - These findings indicate that both inter- 
cluster magnetic interactions and magneto-elastic cou- 
plings are present in general in a specific material, and 
that it only depends on the relative strength whether the 
one or the other is considered as the driving force for the 
phase transitions at the LCs. Hence, as a conclusion, the 
phenomena of the field-induced magneto-elastic instabil- 
ities and field-induced magnetic ordering appear as the 
two extreme sides of one and the same medal — and it 
might be hard to distinguish them from each other in a 
specific material. 



II. BASICS 

In this section, some general results for an effective 
two-level Hamiltonian (TLH) are shortly reviewed. In 
matrix form the TLH is written as 



Htl = 



si A/2\ 
A/2 e 2 ) 



(2) 



where E\ and £ 2 are the bare energies of the two involved 
levels, and A describes a mixing between these levels. 
The approximation of the TLH works well then both the 
gap between the two levels and the temperature are much 
smaller than the gaps to the next higher lying levels in the 
energy spectrum. The energies E± of the two eigenstates 
of the TLH are 



1 
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E± = - (e 2 + £1) ± -v / (e 2 -£i) 2 + A 2 



E 1 ±E a , (3) 



where E + (E_) denotes the energetically upper (lower) 
branch of the avoided LC. The partition function 
Z = 2 exp(— (3Ei) cosh(/3E 2 ) and the free energy F — 
— fcsTln(Z) are functions of the temperature T, and fur- 
ther variables, which are abbreviated symbolically by x, 
i.e., Z = Z(T,x) and F = F(T,x). k B is the Boltzmann 
constant and (3 = 1/ {ksT). For some relevant derivatives 
one obtains 



dF{T,x) 
dx 



d 2 F(T, : 
dT 2 



8E 1 , mj? ,dE 2 
— -tanh(/3£ 2 ) — 

— + [l-tanh(/3£ 2 )]— ^, (4a) 
-4/3 3 [l-tanh 2 (/3£; 2 )](£; 2 ) 2 . (4b) 



A formulation of the expressions in terms of a factor 
[1 — tanh™(/3i? 2 )] (with integer n) is convenient for a dis- 
cussion of the temperature dependence: For T — > one 
finds [1 — tanh"(/3i? 2 )] — > 0, and in the high-temperature 
limit T -> 00 one finds [1 - tanh n (/3E 2 )} 1. E 2 is half 
of the energy gap between the two levels, E + — E- = 2E 2 . 



III. TWO-LEVEL HAMILTONIAN 

In this section, the effective TLH, which describes the 
behavior of a (rigid) spin cluster at a field-induced LC, 
is derived in first-order perturbation theory. The spin 
Hamiltonian H of the system is given by Eq. |T]). The 
type of LCs considered in this work emerge from dom- 
inant Heisenberg interactions and the Zeeman term in 
the cluster, such that the effects of the term Hi can be 
treated perturbationally. For simplicity it is further as- 
sumed that the cluster anisotropy is uniaxial. Then the 
magnetic field is completely specified by its magnitude 
B and its angle (f with respect to the uniaxial z axis. 
No further restrictions are imposed on the treatment in 
this and the next section. It is noted that for most AFM 
wheels, and CsFes in particular, the anisotropy is very 
well characterized as uniaxial^^^ii 4 ^^ 

The most important term in Hi is usually a single-ion 
anisotropy term, 



H 



D 



= J2 Si • • Si 



but also a Dzyaloshinsky-Moriya term, 



H. 



DM 



5> 



[Si x Sj), 



(5) 



(6) 



or a spin-phonon coupling term, etc., may be included, 
depending on the situation under consideration. Since 
the single-ion anisotropy is of general importance, it is 
convenient to introduce a H[ via Hi = Ho + H[. In the 
following and Sec. IV, H[ is mostly disregarded, but it 
will become crucial in Sees. V and VI. 

Since the main part of the microscopic Hamiltonian H 
is isotropic, the appropriate zero-order states are eigen- 
states of the total spin operator S, and classified by the 
quantum numbers S and M. The two states involved in 
the LCs have spin quantum numbers S differing by one 
unit and M = — S*.— The wavefunctions may be writ- 
ten as I7, S, —S) and I7', S + 1, — S— 1), where 7 denotes 
additional quantum numbers. The obvious shorthand 
notation \S) and \S + 1) will be used in the following. 
In first-order the TLH Eq. @ is obtained with Ei = £5, 
£2 = £5+1, and 



£ S = (S\H\S), 
A/2 = (S\H\S+l). 



(7a) 
(7b) 



It is important to note that the quantization axis for the 
magnetic quantum number M is along the direction of 
the magnetic field, and not for instance along the uni- 
axial z axisi 6 ' 27 i 49 ' 55 i 56 For the calculation of the matrix 
elements, an operator such as Sf z , which is expressed in 
a reference frame with the quantization axis along the z 
axis, thus has first to be rotated to a frame with the z' 
axis (the "new" z axis) along the magnetic field. 

For the diagonal matrix elements one directly calcu- 
lates e s = A° s -bS+(S\Hi \S) where A^ is the energy 
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of the level \ S) in zero magnetic field and without H%, and 
the abbreviation b = gu B B is introduced. The first-order 
approximation thus results in a linear field-dependence of 
the bare levels near a LC. The non-diagonal matrix ele- 
ment reduces to A/2 = (S\H\\S + 1), since the Heisen- 
berg and Zeman terms do not mix states with different 
S. As an important result it follows that, in first-order, 
A does not depend on the magnetic field, but only on 
the angle <p, i.e., A = A(<p). 

Because of the importance of the single-ion anisotropy, 
the situation with Hi = Hu is discussed in some 
detail. The matrix elements of Ho were calculated 
previously^ ^ 24 ' 27 ! 49 ' 52 For the diagonal elements holds 
(S\H D \S) = D S S(S - l/2)(cos 2 (^ - 1/3), where D s is 
the zero-field-splitting parameter of the spin multiplct to 
which \S) belongs. One thus obtains 



e s (b,(p) = A s (p)-bS, 



(8) 



where As ((f) denotes the zero-field energy of the level \S) 
due to the combined effects of the Heisenberg interactions 
and the single-ion anisotropy Hu ■ ft is interesting to note 
that the dependence on the angle tp comes entirely from 
the zero- field gaps As(f). 

The magnetic field bo, at which the LC occurs, is 
determined by the condition £i(&o) — £2(^0)1 yielding 
b (<p) = A 2 (<p) - Ai((p) or 



bo(ip) — a + 6(cos 2 <p — 1/3), 



(9) 



with the constants a — A° — A? and b — Ds+i(S + 
1)(S + 1/2) - D S S(S - l/2). 6 -^ The LC field b in gen- 
eral exhibits an angle dependence due to the anisotropic 
terms in Hi. For Ho one finds the generic behavior 
bo(ip) (x (cos 2 ip — 1/3) + const. This angle dependence 
of the LC fields has been observed in high-field torque 
experiments on several molecular AFM wheels, and was 
used to experimentally determine the zero-field-splitting 
parameters £) s ,SiL2i2iI£ 

Concerning the non-diagonal elements of Hp, one finds 

(S\H D \S+l) = cos(psm(p^2 Di(S\f[ 2) {i)\S + 1), 



where Di is the projection of the single-ion anisotropy 
T>i on the uniaxial axis, and Tf (i) the irreducible ten- 
sor operator (ITO) of rank 2 related to the ith spin 
centerj 50 i 51 ' 52 Accordingly, the angle dependence of the 
mixing parameter follows 



A(tp) cx cos ip sin ip, 



(10) 



so that the level mixing induced by Ho is zero for parallel 
and perpendicular fields, <p = 0° and 90°, respectively, 
and maximal for ip = 45°. 

The above first-order results are often a very good ap- 
proximation. In systems, however, with a large mag- 
netic anisotropy, such as the CsFes wheel^^i the effects 
of higher-order contributions become non-negligible and 



would have to be considered for a fully quantitative de- 
scription. For instance, in second-order an additional 
quadratic field term oc b 2 would arise. But even then, 
the first-order results are useful for qualitative or semi- 
quantitative considerations. 

The inclusion of further magnetic anisotropic terms 
via H'i produces additional contributions to E\, £2, and 
A, i.e., adds the respective matrix elements of H[ to 
Eqs. (JSj) and (fTO]) . One term of potential importance is 
the Dzyaloshinsky-Moriya interaction Hom ■ 24,42 For this 
term the diagonal contributions are zero, so that e\ and 
£2 [as well as bo(ip)] are not affected. For the non-diagonal 
element one finds 

(S\H DM \S + 1) = / cos ipdij x + djjy 

-ism<pd ijz )(S\fW(ij)\S+l), 

where f^\ij) is an ITO of rank 132L22 The contri- 
bution of H]jm to the angle dependence of the mixing 
parameter A thus is quite complicated in general. For a 
planar uniaxial cluster with sufficiently high symmetry, 
such as molecular wheels or grids, only dij Z is non-zero 
and A will have a component which varies as sin^i^ 



IV. THERMODYNAMICS 

The thermodynamic properties for magnetic fields 
close to a LC can be calculated analytically, keeping 
in mind Eq. ([8]) and A = A(ip). Since a system with 
uniaxial magnetic anisotropy is considered, the free en- 
ergy F is a function of the temperature T, of the mag- 
nitude of the magnetic field B (or b), and of its ori- 
entation tp: F = F(T,b,(p). In this section, the mag- 
netic specific heat, the magnetization, and the magnetic 
torque is calculated. For the calculation it is noted that 
£2 ~ £1 = bo(ip) ~~ b. It is also convenient to introduce the 
two functions 



G(b,<p) 
H(b,<p) 



b-b ((p) 



2V[6-6oM] 2 + A 2 (^)' 
2V[6-M^)] 2 + A 2 (^)' 



(11a) 
(lib) 



As function of field, G(b, <p) describes a step from -1/2 to 
1/2 centered at the LC field bo, which is broadened by 
the level mixing due to A, as shown in Fig. QJa). The 
field derivative dG/db exhibits a peak at the LC with 
a height of 1/(2|A|) and a full-width-at-half-maximum 
(FWHM) of 2V2 2 / 3 - 1|A| = 1.124|A|. 48 H(b,ip) is es- 
sentially the square-root of a Lorentz function and ac- 
cordingly describes a peak at the LC field bo with a height 
of ±1/2, depending on the sign of A, and a FWHM 
of 2V^|A| = 3.464|A|^ see Fig. [lib). Both G and 
H can be expressed in terms of the reduced variable 
x = {b - bo)/\A\: G(x) = xH(x), H(x) = (2Vx T TT)- 1 , 
which makes the dependencies on 60 and A transparent. 
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FIG. 1: (Color online) Field dependence of the functions (a) 
G(b, tp) and (b) H(b, tp) discussed in the text for various values 
of the level-mixing parameter A. 



The magnetic specific heat is given 
-Td 2 F/dT 2 , which with Eq. (gb| yields 



by cm 



cm (b, tp) 



El 



k B T 2 cosli '(j3E 2 ) 

2 



exp(/32£ 2 ) 



2E 2 



k B T J [l + exp(/32£; 2 )] s 



(12) 



Since 2E 2 is the energy gap between the two levels, the 
specific heat exhibits a Schottky anomaly as function 
of temperature (with a maximum of height 0.4392fcs at 
k B T = 0.4168|2i5 2 |), and a double-peak feature as func- 
tion of field, as was noted and discussed before by several 
authors, to whose works we hence refer to . 23 ' 56 i 58 ' 59 i 60 

The magnetization is given by m = —DF/db (in units 
of gnu) and calculated to 



m(b, ip) = S + - + G{b, tp) tanh(/3£b). 



(13) 



It is convenient to separate the magnetization as m = 
+ 5m, where m<g = S is the zero-temperature mag- 
netization at fields well below the LC field, b <C bo. 5m 
is the change in the magnetization upon passing the LC 
coming from fields below 6q : 



5m(b, tp) 



: + G(b, tp) tanh(/3£ 2 ) 

\ + G{b,<p) - [l-tanh(/?£ 2 )]G(&,<p). 



8m(b) 
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FIG. 2: (Color online) Field dependence of the magnetization 
change Sm(b, tp) for various temperatures. The level-mixing 
parameter was set to A = 0.25. 



(14) 

As function of field, the magnetization exhibits a step at 
the LC, where it changes from S to S +1 (or by one unit 
<7/Lts). Figure [2] shows 5m(b) at various temperatures for 
A = 0.25. The magnetization can be expressed in terms 
of the reduced variables x = (6— 6 )/|A| and t = k B T/\A\ 
[(3E 2 = V x 2 + 1/ (2t)]. In these units, the variation of the 
behavior with 6 an d A becomes most apparent. 

At low temperatures, k B T -C |A|, the magnetization 
change is given by 1/2 + G, and is independent on tem- 
perature. At high temperatures, k B T ^> |A|, it reduces 
to 1/2 + tanh[(6 - b Q )/(2k B T)}/2 as for a true LC, i.e., 
the effects of level mixing become indiscernible in this 
temperature regime. 

The magnetization step is broadened by the effects of 
temperature and non-zero level mixing A. The broad- 
ening may be analyzed by considering the width of the 
peak in the field derivative dm/ db. At low temperatures, 
k B T <C |A|, the peak is characterized by a temperature- 
independent FWHM of 1.124|A|, and at high tempera- 
tures, k B T 3> |A|, by a linearly increasing FWHM of 
41n(\/2 + l)k B T = 3.525fc B T.£2iii£ With temperature, 
also the shape of the peak changes, so that a careful anal- 
ysis of the shape of experimental dm/ dB curves in prin- 
ciple provides additional information on the level mixing. 

The magnetic torque is given by r = dF/dtp^ and 
calculated tc 



J>2 



r(b, tp) 



1 / gAi 

2 V dtp 



dA 2 



dtp 



tanh(^£; 2 
<9A 



(15) 



This equation can be transformed into a more useful form 
by considering i) (A x +A 2 )/2 = Ai +b /2 and Eq. fli]). 
and by separating ii) the torque into t = T<g + 5t, where 
T<g = dAi/dtp is the torque for fields well below the LC 
field, and St the change in the torque upon passing the 
LC. One finds^ 



5r(b, tp) = Sm(b, tp) 



^ + ff(6,^)^tanh(/3£; 2 ) 
dtp dtp 
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-{l-t&nh(f3E 2 )}H(b, tp)^. 



(16) 



Interestingly, the torque consists of two contributions, 
5t = 5ti + St2, where 8t\ refers to the first term on the 
r.h.s. of Eq. (Til))) , and St^ to the remaining terms. The 
first contribution, 8t±, is proportional to the magnetiza- 
tion change 5m. It accordingly also produces a step as 
function of field in the torque signal, which is broadened 
by a non-zero level mixing and temperature. The second 
contribution, St2, as it is proportional to H(b,tp), pro- 
duces a peak in the torque signal centered at the LC field. 
The field dependence of St 2 is depicted in Fig. [3] for dif- 
ferent temperatures and several values of A. At low tem- 
peratures, ksT <c |A|, the peak follows H(b,tp) and is 
accordingly characterized by a temperature-independent 
height of (l/2)dA/dtp and FWHM of 3.646|A|4£ With 
increasing temperature the peak smears out to disappear 
at temperatures ksT ^> |A|. As with the magnetization, 
the torque is expressible in terms of the reduced vari- 
ables x and t, clarifying the trends with A and T visible 
in Fig. [31 Several points shall be noted: 

1) For the special case of zero level mixing (A = 0), 
most of the results have been obtained before, in par- 
ticular that both the field-dependent magnetization and 
torque curves exhibit steps at the LC field. In fact, from 
Eq. (|16[) one finds 5r(b, p) cx 8m(b, tp), since 5r2(b, tp) = 
for A = 0. In reality, however, A never will be exactly 
zero, but as long as it is sufficiently small, the torque 
will be characterized by steps [small means OA/ dtp <C 
dbo/dtp, which is obeyed if A(tp) is either small in mag- 
nitude or almost constant as function of angle]. Vice 
versa, whenever r(6, tp) is observed to be mostly step- 
like, one can conclude that (i) A is very small and (ii) 
5r(b, ip) oc Sm(b, tp). This situation has been observed for 
instance in many AFM w heels4iSiSiLSi2ii2iSi5£ 

2) If level mixing, however, is significant, then the mag- 
netization and torque curves in general will be very dif- 
ferent, since then the peak-like contribution <5t2(6, tp) sig- 
nificantly adds to the torque signal. A similar contribu- 
tion is not present in the magnetization, and it is thus a 
unique feature of torque. This contribution is rooted in 
the fact that the mixing parameter A does not depend 
on the magnetic field but on the angle tp, and is hence 
directly probed by torque. If level mixing is actually very 
strong, i.e., dA/dtp S> dbo/dtp, then the torque curve will 
be characterized by peaks at the LCs. At nonzero tem- 
perature and/or strong level mixing the peaks of several 
nearby LCs may superimpose to yield an oscillatory field 
dependence of the torque^ This situation has been ob- 
served for the Mn-[3 x 3] grid and the modified wheels 
Cr 7 Ni and Cr 7 Zn^>^ 

3) Interestingly, if one writes the ground-state wave- 
function as 
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FIG. 3: (Color online) Field dependence of the torque contri- 
bution 5t2 (b, tp) at several temperatures for the level-mixing 
parameters (a) A = 0.1, (b) A = 0.25, and (c) A = 0.5. 



with 



,s 2 = 1, one finds 



s 

cs 



l/2 + G(b,tp), 
H(b,tp). 



(18a) 
(18b) 



|0) =c\S)+a\S + l), 



(17) 



Hence, at zero temperature, 5m(b,tp) and <5ti(6, p) di- 
rectly reflect s 2 , while Sr2(b,tp) reflects the product cs. 

4) From the above points it is clear that the peaks in 
the torque curves due to 5r2(b,tp) are a direct signature 
of level mixing, as pointed out previously in Ref. [26l In 
this work it has been also shown that the peaks in the 
torque are related to the enhancement of the quantum 
fluctuations of the total spin S by a level mixing. In fact, 
a derivation has been given which yielded r oc 2BAS Z , 
where (AS^) 2 = (5 2 ) — (S z ) 2 denotes the quantum fluc- 
tuations in S z .— This relation cannot be correct, or com- 
plete, since for zero level mixing it would predict r = 0, 
in contrast to the theoretically expected and experimen- 
tally observed step-like torque curves. In view of the 
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above findings, however, 2BAS Z can be associated with 
the contribution <5r 2 , as both are related to cs. 

5) The generic behavior of the angle dependence of 
the LC field is bo ((p) oc (cos 2 <p — 1/3) + const, so that 
dba/dtp oc sin ip cos p. Thus, the step-like part of the 
torque will follow 5Ti(b,tp) oc Sm(b) sin tpcos ip, which is 
the expected angle dependence. 

6) In the above treatment an isotropic g factor was 
assumed. The calculations for anisotropic g factors are 
simple but the results are lengthy. As a general trend, an 
anisotropic g factor produces an additional contribution 
to the torque (and only to the torque) which increases 
linearly with the magnetic field in first approximation. 53 

V. MAGNETO-ELASTIC HAMILTONIAN 

Within the spin-Hamiltonian formalism a magneto- 
elastic (ME) coupling is introduced by allowing the 
magnetic parameters (coupling constants, anisotropy pa- 
rameters, g factors, etc.) to depend on the distor- 
tion coordinates Q k = — u^o of the atoms in the 
lattice^iiSiiSiHieSieeie! u ^ q ig t ^ e equil i br j um and Ufc 

the distorted position vector of the kth atom; the Q fc 
are usually normal coordinates. For the exchange cou- 
pling constants, for instance, hold Jy (Q 1; Q 2 , . . .) = 
Jij + S Jij(Qx, Q2, • • •)■ Assuming small deviations, a 
Taylor expansion yields 5Jij (Q 1; Q 2 , . . .) ~ J2k a ij,k ' u k, 
with the ME coupling constants o^fe = dJij/duk- In 
the following, it is assumed that only one distortion 
mode, and hence only one (scalar) distortion coordinate 
Q, is relevant. For the exchange constant one then finds 
JijiQ) = Jij + 5Jij(Q) with S Jij — ctijQ. Similar rela- 
tions hold for the other magnetic parameters D^, d,j, 
and so on. 

Inserting the modulated magnetic parameters into the 
microscopic spin Hamiltonian produces additional, Q- 
dependent terms which will be denoted collectively as 
Hq. For the terms considered in this work, one finds 

Hq = - ) SJjjSj ■ Sj + Si • &Dj • Si 

i=£j i 

+ ^<Sd, r (S,xS J ) + ..., (19) 

The combined system of the magnetic molecule and 
the elastic lattice is then described by the spin-phonon 
Hamiltonian 

H sp ^ ph =H + f+^Q 2 , (20) 

where Hq is included in H via Hi or H[, respectively. 
T is the kinetic energy of the phonons and k the elastic 
constant. 

In this work, the phonons are treated as classical oscil- 
lators, which is justified if the lattice dynamics is much 
slower than the spin dynamics. This situation is known 
as the adiabatic limit and has been often used to infer the 



ground state in spin-phonon systems i ' i 65 Under this 
conditions, the ground state is obtained for zero kinetic 
energy of the phonons (atoms at rest) and by minimizing 
the potential 

V(Q) = E(Q) + ^Q 2 (21) 

of the total system, where E(Q) is the ground-state en- 
ergy of the magnetic part H (which now depends on 
Q due to the inclusion of Hq). One may approxi- 
mate E(Q) = —aQ K , with some positive constants a 
and k. A stable minimum is obtained for k < 2 at 
Qo = 2 ~^J na/k > 0. This situation corresponds to un- 
conditional instability. For k = 2 conditional instability 
is realized, where V(Q) < is obtained for any Q > if 
a > k/2. For k > 2 the only solution is Q = 0, i.e., no 
distortion is obtained. 

Within the framework of the TLH approximation, 
E(Q) = E_(Q) (_E_ depends on Q via £1, £ 2 , and 
A). Apparently, the ME coupling enters in two ways, 
namely via the diagonal matrix elements (S\Hq\S) , and 
via the non-diagonal matrix element (S\Hq\S + 1). The 
first mechanism shall be called diagonal, the second non- 
diagonal ME coupling. These two coupling modes have 
to be carefully distinguished. 

Before proceeding further, a point concerning the 
Heisenberg-exchange part shall be clarified. Within 
the context of the present work it sets the largest 
energy scale. In several models treated in the lit- 
erature it is actually the only term in the spin 
Hamiltonianj 41 i 54 ' 68 i 69 ' 70 i 71 Two cases have to be distin- 
guished, namely those where the Heisenberg exchange 
results in (i) a degenerate or (ii) non-degenerate ground 
state. Case (i) is encountered, for instance, in spin trian- 
gles or tetrahedrai 68 i 71 The description bears some anal- 
ogy to the Jahn- Teller effect ( JTE) , and the phenomenon 
hence has been denoted frequently as the spin- JTE. Such 
situations have been analyzed theoretically several times 
in the literature ) 68 i 69 ' 70 i 71 but to the best of our knowl- 
edge no experimental evidence has been reported to date. 
Case (ii) is clearly distinct, and is the one of interest 
in the present work. Here, an unconditional lattice in- 
stability is obtained in zero field for gapless spin sys- 
tems, such as a spin-1/2 chain (spin-Peierls effect)^ 4 - For 
gapped systems, such as the finite spin clusters consid- 
ered here, a spontaneous instability typically does not 
occur. With appropriate magnetic fields a degenerate 
ground state can be created via the LC mechanism, but 
for a purely isotropic model a spontaneous instability 
nevertheless does not occur: The degeneracy at the LC 
points cannot be lifted by whatever distortion since the 
two involved levels are eigenstates of S with different val- 
ues of S. They hence cannot be mixed by the Heisenberg 
exchange. This already indicates the crucial role of a 
magnetic anisotropy in Hi for a spontaneous instability 
at a LC in case-(ii) systems. 

The analysis of the diagonal ME coupling requires the 
evaluation of (S\Hq\S). There is no general reason why 
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(S\Hq\S) should be zero, but a non-zero value would 
lead to somewhat unphysical situations. The effect of a 
diagonal ME coupling would be a shift of the LC fields 
by a distortion, bo = bo(Q)- A change of bo, however, 
corresponds to a change of the zero-field splitting, see 
Eq. JHJ). For the exchange coupling, e.g., this would im- 
ply a change of the average coupling strength, i.e., of 
the overall magnetic energy scale. This is not what one 
normally expects. For instance, for a ring-like model 
the generic distortion mode would be a dimerization, 
upon which the average J2i is n °t altered. Hence 

SJi^i + i — 0, implying a zero diagonal element. Simi- 
lar arguments can be put forward for the other magnetic 
terms (Hdm gives zero anyhow). Thus, for typical physi- 
cal models the diagonal ME coupling is zero in first-order. 
Second-order terms, however, are always expected to be 
non-zero, and one accordingly concludes that 



(S\H Q \S)<xQ 2 



(22) 



corresponding to conditional instability. These consider- 
ations are in agreement with the findings of Refs. 4lll54T |. 
For a gapped quantum spin system one expects the 
generic behavior E(Q) cx Q 2 . 

The situation is different for the non-diagonal ME cou- 
pling (S\Hq\S + 1). For the Heisenberg contribution 
in Hq the matrix element is strictly zero, but for the 
anisotropic terms it will be non-zero in general. This 
is demonstrated by the simplified example of a regu- 
lar AFM spin ring, for which a Dzyaloshinsky-Moriya 
interaction emerges upon a dimerization, dj,i+i, z (Q) = 



(-IfaQ (d 



2,2+1, IE 



= d 



0). Then (S\H Q \S 



1) « Ei(-l)*aQ(S|?i (M + l)\S + 1) oc Q, 



(S\f^(i,i + 1)\S + 1) = {-If (S\T^(1)\S + 1). The 
point is, in a qualitative language, that the diagonal el- 
ements essentially probe the average of the modulation 
(in our example J2i <Wm+M = 0)> while the non-diagonal 
matrix element is sensitive to the span of the modulation 
[in our example l)*^M+i,* ^ 0]- Hence, in general, 



Mi 



(S\H Q \S + 1) cxQ + const, 



(23) 



which allows for unconditional instability. 

The above considerations imply the following TLH to 
describe the ME effects at a LC: 

fr _( e s (b,<p) A(^,Q)/2\ k 2 
HsF -\ Afo,Q)/2 e g+l Q>,<p) 



A(<p,Q)=A(<p,0) + a(<p)Q. 



(24) 



The essential elements of this model are that (i) the diag- 
onal elements do not depend on the distortion Q, (ii) the 
non-diagonal elements do not depend on the magnetic 
field b, and (iii) A(Q) varies linearly with Q. 

The model explicitly allows for the possibility of a 
non-zero non-diagonal matrix element even for the non- 
distorted molecule, i.e, for A(0) ^ 0iZ2 It is easy to show 



that the system exhibits an unconditional instability for 
any value of A(0). However, the most interesting situa- 
tion arises for sufficiently small A(0) (what small means 
will be made precise later) , and this case is also the one of 
relevance for CsFes^ The behavior for small A(0) is not 
essentially different from that for zero A(0), and in the 
following A(0) = is hence always assumed. The case of 
A(0) ^ will be considered at the end of the next sec- 
tion. The other approximations and limitations leading 
to Hsp have been carefully discussed in the above. 

Formally, Hsp is equivalent to the standard TLH dis- 
cussed in the JTE. Indeed, the problems are in many 
respects similar, and much of the insights and results de- 
veloped for the JTE can be directly carried over to the 
current problem (key words to mention are cooperative- 
JTE, pseudo-JTE, etc.)i24 The physical difference, how- 
ever, is that the JTE is related to a degeneracy in the 
electronic system, while here it originates from a degen- 
eracy in the spin system. The effect described by H$p 
hence might be called a spin- JTE. This notation is not 
quite satisfying, since the underlying mechanism is differ- 
ent from the spin-JTE for Heisenberg spin systems with a 
degenerate spin ground state, which was mentioned be- 
fore. The effect discussed here therefore is tentatively 
called field-induced spin-JTE (FISJTE). 

An interesting feature, as compared to the electronic 
and spin-JTE, is the possibility to continuously tune 
through the instability via an applied magnetic field. At 
the LC, es+i — £s is zero, but can be adjusted to any 
value by moving away from the LC point, realizing a sit- 
uation similar to the pseudo-JTE. In the next section, 
the behavior of the system as function of magnetic field 
is considered in detail. 



VI. FIELD-INDUCED MAGNETO-ELASTIC 
INSTABILITY 



For the discussion of the field-induced ME instabil- 
ity it is convenient to first describe the potential energy 
surfaces (PES) produced by Hsp- Two PES are ob- 
tained, V(b, Q)± = E±(b, Q) + (k/2)Q 2 (where of course 
V- = V). The dependence on Q is shown in Fig. Ufa) 
for several magnetic fields. At the LC, b — bo = 0, the 
PES are degenerate at Q — and exhibit two minima at 
Qmax _ a /(2k). With increasing distance from the LC 
field, i.e., increasing \ b— bo\, a gap opens between the PES 
due to the Zeeman splitting, and the two minima accord- 
ingly shift to lower Q values. For fields below bo — b c or 
above bo + b c , where b c = a 2 /(2k), the minima disappear 
completely (in the context of the electronic JTE, b c is the 
JT energy). 

The spontaneous distortion in the ground state is de- 
termined by the minima in the PES, dV/dQ — 0. It 
is more convenient, however, to discuss V(Q) — V(Q) — 
E-(Q = 0), which directly measures the gain in potential 
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FIG. 4: (Color online) (a) Distortion dependence of the two 
PES V±(b, Q) for several fields b — bo and (b) field and distor- 
tion dependence of the potential V(b, Q) — V- (b, Q)—E- (b, 0) 
for positive values of Q (b c = A^ ax = 1, k = 1/2). 



energy due to a distortion. V is determined to 



V(b,Q) 



1 



1 



-V(b-b ) 2 + A(Qr + -\b-b \ + -Q 2 . 



(25) 

The field and distortion dependence of V is shown in 
Fig. HJb) for positive values of Q (and k = 1/2). One 
clearly observes a field range around the LC where V 
is negative for non-zero Q. At these fields the system 
gains energy by distorting to a finite distortion Q, i.e., a 
spontaneous instability occurs. 

From minimizing V(Q), or V(Q), the condition for 
the equilibrium distortion Qo, or the equilibrium mixing 
Ao = A(Qo) = a Qo, respectively, is obtained as 



(b - b ) 2 + A 2 a 



(26) 



This is of the form x 2 + y 2 = r 2 , hence the mixing Ao 
(or Qoi either one can be used as the order parameter) 
describes a semi-circle with radius b c as function of the 
relative magnetic field b — bo, a s shown in Fig. OJa). For 
fields below 6 — b c and above bo + b c , one finds A = 
and the system is undistorted. For fields in the range 
[bo — b c ,bo + b c ] , however, Ao becomes non-zero signaling 
the spontaneous instability. At exactly the LC, the mix- 
ing and the distortion are maximal, assuming the values 
and (AjJ 10 * = aQ nax = b c = a 2 /(2k)). 
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FIG. 5: Field dependence of (a) the equilibrium mixing pa- 
rameter (or order parameter) Ao and (b) the energies E±—Ei 
of the two levels (6 C = A^ ax = 1, k = 1/2). The dashed lines 
in (b) show the energies of an ordinary avoided LC with a 
mixing parameter A = Ao" 1 *. 



Concerning the ground-state wavefunction |0), 
Eq. (jTTJ) , the coefficient s = (5+110) is calculated 
to s 2 = [1 + (b - b )/b c ]/2 by inserting Eq. (|2BJ> into 
Eq. (|18a|) . Thus, in the field range [bo — b c ,bo + b c ] the 
parameter s 2 increases linearly from to 1 as function 
of field (and is and 1, respectively, for fields outside 
this range). 

With Eqs. ^ and (|26p. the energies of the two lev- 
els are calculated to E± = (e 2 + e x ± A nax )/2 for fields 
in [bo — b c , bo + b c ]. Outside this field range, the ener- 
gies exactly correspond to that of a true LC. The field- 
dependence of E± — Ei is shown in Fig. [5fb). Thus, 
the distortion is such as to maintain a minimum gap 
of value Aq™ 021 between the energy levels. For compar- 
ison also the energy diagram of a conventional avoided 
LC with a mixing A™ 01 ' is shown, in order to emphasize 
the distinctly different behavior of the present model. In 
an avoided LC, the mixing is field independent. In the 
present model, in contrast, the mixing is zero for fields 
outside the range [bo — b c , bo + b c ] , and non-zero within it 
due to the instability, see Fig. [H^a) (it is recalled that the 
mixing is directly related to the spontaneous distortion 
via A = aQo)- 

The model permits an analytical calculation of the 
thermodynamic functions at zero temperature (or suf- 
ficiently low temperatures, fceT Ao lax ). The specific 
heat is zero, since the energy levels never come closer 
than A r nax . Concerning the calculation of the magneti- 
zation and torque, the two functions G and H, Eq. (jlip . 
simplify to G = {b - b )/{2b c ) and H = A /(2A nax ) for 
fields in the range [b —b c , b +b c ]. Outside this range they 
become G = ±1/2 and H = 0. Inserting this in Eqs. (CO 
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FIG. 6: Field dependence of (a) the magnetization change 
dm, (b) the field-derivative of 5m, (c) the torque contribution 
St2, and (d) the field derivative of the torque contribution 5t2 



(be 



1, fc = l/2). 



and (|16p (and considering T — > 0), one obtains^ 



5m(b, ip) 
5r(b, ip) 



b-b 



b c 
b-b 



db 1 A dA 
dp 2 A^ ax dip ' 



(27a) 



(27b) 



As function of field, the magnetization change 5m in- 
creases linearly from 5m = at b = bo — b c to 5m = 1 at 
b = bo + b c , as shown in Fig. EJa). Figure [6jb) presents 
the field derivative of the magnetization, dm/ db. 

As emphasized in Sec. IV, the torque consists of two 
contributions, 5t = 5t% + 5t2- The first contribution is 
simply proportional to the magnetization change, 5ti = 
5m(db / dip) , see Eq. ((Tfj)) . and hence also increases lin- 
early in the field range [&o — b c , bo + b c ] . The field de- 
pendence of the second contribution needs more consid- 
eration. As will be discussed in a moment, the term 
dAo/dp should be interpreted as (da/dip)Qo- With re- 
spect to the field dependence it is hence proportional to 
Ao. This yields 5t2 <x Aq. The field dependence of 5ti 
and its field derivative are shown in Figs.^c) and[5Jd), 
respectively. 5t2 displays a dome-shaped behavior in the 
field range [&o — b c , bo + b c ]. Both types of torque curves, 
"linear slope" and "dome shaped", were observed in ex- 
periments on the AFM wheel CsFes^ 

In general, since 5t = 5t\ + 5ti , the torque is the sum 
of a linear slope and a dome-shaped contribution. The 
respective weights are governed by the factors dbo/dip 
and (da/ dip) /a, or the ratio of them, respectively. The 
actually observed torque profiles hence will depend on 
the angle ip in general, as well as on the particular sys- 
tem under consideration. A careful investigation of the 



angle dependence of the torque profiles thus can provide 
information on the relevant ME coupling terms in Hq. 

At this point an apparent contradiction shall be ad- 
dressed: In the field range [bo — b c , bo + b c ] the energies 
E± as well as the magnetization exhibit a linear slope, 
in conflict with the relation m = —dE/db valid at zero 
temperature. Also, the above procedure of just inserting 
Eq. ([2S]) into the equations of Sec. IV seems question- 
able, since A is implicitly field (and angle) dependent 
via Eq. (f2"rj)) . This violates the assumption dA/db = in 
the derivations of Sec. IV. However, strictly spoken the 
PES V± should be used in the calculations, and not E±, 
since the combined system of spins and lattice vibrations 
is considered. This would add terms due to dQ^/db and 
dQo/dip to the magnetization and torque. Such a more 
rigorous calculation, however, yields exactly the same re- 
sults as with using the equations of Sec. IV and neglect- 
ing the implicit field and angle dependence of Ao. That 
is, for the thermodynamic functions Ao = aQo can be 
used in the sense of dQo/db = dQo/dp = 0. The en- 
ergy spectrum shown in Fig. Efb) is not suited to derive 
the thermodynamics, but correctly reflects the gap as it 
would be observed in spectroscopic measurements. 

The behavior for non-zero temperatures may be 
treated along the lines of, for instance, Ref. [69|. With 
the two PES V±, the free energy F = F(T,b,ip,Q) be- 
comes 



F = Ei 



In 







2 cosh T;V(b-b ) 2 + (aQy 



(28) 

The equilibrium distortion Qo is determined by dF/dQ = 
0. It is convenient to switch to the order parameter A = 
aQ, and to express the resulting condition in terms of the 
reduced variables b = (b — bo)/b c , Ao = Ao/A r nax , and 
T = T/T™ ax (where k B T" lax = b c /2; we recall A^ = 
b r ). In these units, one obtains 



tanh(V/T) = X, 

x^Jbi + Al 



(29a) 
(29b) 



This sort of transcendent equation often appear in mean- 
field theories, for instance the Weiss molecular-field the- 
ory of ferromagnetsi^ This shows that the above theory 
actually establishes the mean-field theory of the FISJTE. 

At the LC, where 6 = and hence X oc Aq, the graph- 
ical solution of Eq. (|2"9"|) proceeds exactly as in textbooks 
on the Weiss ferromagneti^ The critical temperature, 
below which a spontaneous distortion occurs, is deter- 
mined by tanh(V/r c ) sa X/T c = X, which yields f c = 1 
or k B T c = k B T™ ax = bc/2, respectively. At the LC, 
T c assumes its maximal value, T™ ax . For fields below 
or above the LC, T c is reduced, in accord with the find- 
ing that the distortion is largest at the LC point at zero 
temperature [see Fig. [5fa)]. 

The zero-temperature behavior is easily rederived. For 
T = 0, Eq. (JUJ) reduces to the condition b 2 + A^ = 1, 
which is equivalent to Eq. 
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FIG. 7: (Color online) Field dependence of the reduced criti- 
cal temperature T c (solid line) . The dashed (green online) and 
dash-dotted (blue online) lines represent the approximate so- 
lutions b » ^/3(1 - T c ) for |6| < 1 and f c » 2 ln[2/(l - I&I)]" 1 
for b — > ±1, respectively. 
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FIG. 8: (Color online) Field dependence of the order pa- 
rameter Ao for the indicated values of the ratio A(0)/Ao la:c 
(b c = AJ 1 " = 1, k = 1/2). 



The critical temperature for an arbitrary field is given 
by the solution of Eq. (29]) for A = 0, i.e., tanh(6/T c ) = 
b, which yields 



T c (b) = 2b In 



1 + 6 



1 



(30) 



The dependence of T c on b (or 6, respectively) is displayed 
in Fig. [3 Near the LC point, where |6| -C 1, one finds 
b « \/3(l — T c ), corresponding to a critical exponent of 
1/2 as expected for a mean- field theory. Near the upper 
or lower critical fields, b — > ±1, the critical temperature 
f c goes to zero according to T c w 21n[2/(l - |6|)] _1 . The 
field dependencies of these approximations are also dis- 
played in Fig. [7l 

The analysis so far was based on the assumption of 
A(Q) = aQ or A(0) — 0, respectively. In real materials, 
however, A(Q) = A(0) + aQ could be more appropriate. 
For instance, in the AFM wheel CsFcs the molecular C4 
symmetry permits next-neighbor Dzyaloshinsky-Moriya 
interactions, which would result in A(0) 7^ [for the 
specific example of CsFes, however, experiments indicate 
that A(0) is in fact rather small 73 ]. The effect of a non- 



zero value of A(0) at zero temperature shall be shortly 
discussed. 

The minimization of V{Q) as in the above, but with 
A(Q) = A(0) + aQ, reveals that the system is uncondi- 
tionally instable for any magnetic field, and not only in 
the range [60 — fee, 60 + 6 C ] as in the case of A(0) = 0. This 
is expected since the unconditional instability is not re- 
lated to A(0) — 0, but to the presence of a non-diagonal 
ME coupling. An analytic solution is not possible for 
general values of A(0). Figure [8] hence displays the nu- 
merically calculated field dependencies of the order pa- 
rameter Ao for various values of the ratio A(0)/ A™ 01 . 
For A(0) = of course the previous result is reproduced. 
Although the system is unconditionally instable in the 
whole field range for any value |A(0)| > 0, Fig. [5] reveals 
that for small A(0), i.e., A(0) < A^ ax , the distortion 
outside of the field range [60 — be, bo + b c ] is very small as 
compared to the maximum distortion Q™ ax — A Q nax /a. 
The behavior remains very similar to the one obtained for 
A(0) = 0, only the onsets of the distortion at the fields 
bo — b c and bo + b c are less sharp. Hence, the magnetiza- 
tion and torque curves will look similar as in Fig. [SJ but 
with the sharp features at bo — b c and bo + b c rounded by 
a small A(0). 

For large values of A(0), A(0) > A nax , the situa- 
tion changes. In this limit the minimization of V(Q) 
yields for the equilibrium distortion the simple result 
Qo = (oi/k)H(b), where in the function H(b), Eq. (Tib), 
one has to insert A(0) for the level-mixing parameter A. 
In terms of the order parameter Ao this corresponds to 
Ao = 2A Q nax H(b), which shows that Ao will not exceed 
Ag™ 031 , and in particular will remain much smaller than 
A(0). Physically — it is now more intuitive to think in 
terms of a small A™ ax or weak ME coupling, respectively 
— the coupling of the spin system to the lattice is so weak 
that a lattice distortion does not affect the energy of the 
system much. The distortion traces the (field-dependent) 
energy gap, as encoded in H(b), but does not modify it 
significantly. The magnetization and torque curves thus 
essentially show the profiles found in Sec. IV for zero ME 
coupling, with a broadening parameter A = A(0). Ac- 
cordingly, for A(0) ^> A™ 1100 a ME coupling will not lead 
to pronounced effects. This is reasonable, it just means 
that in order to raise the effects of the field-induced ME 
instability to an observable level, the system should ex- 
hibit a sufficiently strong coupling to the lattice and/or 
sufficiently soft vibration mode. 
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